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Abstract 

Several liters of an organic liquid scintillator (LS) loaded with Nd have been 
made. We report on performances of this scintillator in terms of optical 

1 Ph properties, radiopurity and light yield for a Nd concentration of 6.5 g/1. A 

possible application to search for the 150 Nd neutrinoless double beta 
decay with a 10-ton scale LS detector is discussed together with further 
(X) improvements. 
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1. Introduction 

In this paper we report on a number of measurements performed to char- 
^ acterize the properties of an organic liquid scintillator based on pseudoc- 

umene and loaded with Nd. This work has been carried out within the 
MetaLS project pp. In particular, we will report on optical measurements 
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made in order to determine the scintillator response and on radiopurity mea- 
surements and purification methods, to discuss possible applications in the 
double beta decay searches field. 

Several liters of Nd-loaded organic liquid scintillator (NdLS) have been 
produced following the methodic developed to produce Yb, In and Gd loaded 
scintillators [2], [3], [1], 0,10], [I]- Measurements were carried out in small sam- 
ples and in a 2.5 1,1m long rectangular quartz cell equipped with two 3-inch 
photomultipliers. Radiopurity measurements were performed by Inductively 
Coupled Plasma Mass Spectrometry (ICP-MS) and 7-spectroscopy by High 
Purity Germanium detector (HP-Ge) facilities at the Gran Sasso Laboratory. 

We have also studied and report here on the performances of a 10-ton 
scale liquid scintillator detector loaded with Nd to search for the 150 Nd 0^/3/3 
decay. 

2. The Nd-loaded organic liquid scintillator and its optical prop- 
erties 

A 2.5 1 batch of NdLS based on Nd carboxilate (Nd-CBX) having a Nd 
concentration ([Nd]) of 6.5 g/1 has been prepared in spring 2008 starting 
from a concentrated master solution ([Nd] ~ 50 g/1) prepared in December 
2007. The adopted Nd concentration is the result of a trade-off between the 
Nd target mass, the light yield and the attenuation length of the doped LS. 
The choice of the carboxylic acid, the proper reagent purification procedures 
as well as the Nd-CBX solid salt synthesis and drying recipe, are the results 
of a R&D activity carried on in the years 2006-2007 and are a follow-up of 
the know-how developed in our working group to produce large quantities of 
stable and highly performing metal (Gd and In) loaded liquid scintillators. 
The ability in producing highly Nd concentrated master solution allows to 
easily produce NdLS batches of several tens of liters in a chemical research 
laboratory environment. 

The light Yield (LY) of the prepared NdLS and of the reference LS (pseu- 
documene based), both using PPO as primary fluor at concentration of 1.5 
g/1, was measured in small cylindrical vials (5 ml), optically coupled to a 
PMT window by silicon grease. The vials were wrapped by teflon tape to 
improve the light collection. As reference the LY of PC with the same fluor 
concentration, measured in the same vial, is adopted (~ 10000 ph/MeV). The 
vial is irradiated by a 137 Cs source (E 7 = 662 keV). The measured LY values 
come from the comparison of the NdLS to the reference LS Compton spectra 
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as shown in FigjTJ The measured light yield, at [Nd]=6.5 g/1 and [PPO]=1.5 
g/1, is ~ 75% of that of PC based LS at the same fluor concentration. 




Figure 1: LY of NdLS for different Nd 
concentration of PPO and of the ref- 
erence LS as measured in small vials. 
For LS 1.5 g/1 of PPO as primary fluor 
is adopted, while for NdLS BPO at 4.0 
g/1. The use of BPO at the adopted 
concentration leads to a LY larger of 
10-20% than with PPO at 1.5 g/1. 



Figure 2: Emission of NdLS with PPO(3g/l) and 
BPO(3g/l); excitation 267 nm. Nd-CBX in PC 
absorbance spectrum is superimposed. 



The optical characterization of the 2.5 1 NdLS batch is completed by a 
spectrophotometric measurement of the light absorption spectra in the rele- 
vant wavelength range (300-450 nm) and by fluorimetric measurement of the 
light emission spectra. A challenge is the presence of Nd characteristic ab- 
sorption bands, and of other colored rare earths that usually come with Nd, in 
the wavelength region where PMTs quantum efficiency is maximal; this cause 
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scintillation light loss. For this reason, particular attention has been put on 
the purification from optically active and from radioactive contaminants in 
view of the specific scientific application (Ozz/3/3 decay search). Adopting our 
reagent purification procedure, we are able to obtain a transparent region 
around 400 nm as shown in Fig|2j In the same figure, the emission spectra 
of NdLS prepared with the same concentrations of either PPO or BPO as 
primary fluors, are superimposed; the emission spectra match quite well the 
NdLS transparent region. 

3. Light Yield and attenuation length measurement 

Besides the measurements described above we have also performed, start- 
ing in spring 2008, tests with a 1-m long rectangular quartz cell filled with 
scintillator. The cell has been filled first with undoped PC-based (LS) scin- 
tillator (PC, [PPO] = 1.5 g/1) as reference, then with NdLS (Nd-CBX in PC, 
[Nd]=6.5 g/1, [PPO] = 1.5 g/l,[MSB]= 30 mg/1). The goal is to measure on a 
few liters scale: 

• the light yield to be compared to the value obtained from several 
milliliters samples 

• to measure the attenuation length which is a fundamental parameter 
for any application of the NdLS. 

The cell is equipped with two 3-inch photomultipliers (PMTs) optically cou- 
pled to the quartz by silicon grease at the two end sides. The cell is wrapped 
with a VM2000 sheet to increase the light collection. The two PMTs are 
put in coincidence and their charge values are stored. The cell is irradiated, 
scanning it along its main axis, by a collimated 100 kBq 137 Cs source. For 
each source position, we build up the Compton spectrum for each PMT; the 
half-height of the Compton shoulder, L^E, is used to compare data with and 
without the Nd. The scintillator light attenuation length is then estimated, 
plotting the quantity L^E against the relative distance to one PMT and 
fitting with the function: 

y(x) = ^ + k 2 e~^ (1) 

The first term reproduces the direct light collection, relevant close to the 
PMTs, while the second the light collection attenuated by the NdLS. ki^ 
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Figure 3: Charge spectra of both PMTs in the 1-m cell test when the 13 Cs source is 
loacated at the center. 



and A are unknowns and determined from the data. In Fig. [3] we show the 
charge spectra of the two PMTs when the source is located at the center of 
the cell. 

In Fig. [4] and [5] we show the results of measurements performed in 2008 
with and without Nd, respectively. When fitting the data sets with the 
function (JTJ, we determine: 

A e J£ = 168 ± 12 

\eff _ 74 _|_ o 
IV NdLS 

As the Borexino grade PC LS attenuation length is known to be of the 
order of ~8m [8], and we measure it to be ~8m when measuring by spec- 
trophotometry in 10cm long cell prior the lm cell PC filling, it is evident 
that the value obtained from the lm cell measurements is affected by geo- 
metric effects related to light propagation. Therefore to extract the NdLS 
light attenuation length from our data we adopt the following procedure: 

1. we normalize the NdLS L^E to the PC Ljy| and we plot it as a function 
of the PMT to source distance. It is shown in Fig. [6} The exponential 
decay trend can now be attributed to the light absorbance of the NdLS. 

2. we determine the effective cell optical path taking into account the light 
reflection at the external quartz surface: after equalizing in timing and 
in amplitude the PMT signals when the source is placed at the cell 
central position, we then move the source at one end. From PMT 
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Data 2008 - lm cell 
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Figure 4: Distribution of L^E for both 
PMTs against the relative distance to 
PMT2 (see text for detailes) for LS wi- 
htout Nd. 2008 data. 



Data 2008 - lm cell 
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Figure 5: Distribution of Lf,f for both PMTs 
against the relative distance to PMT2 (see text 
for detailes) for LS wiht Nd. 2008 data. 



signal timing measurements we can evaluate the effective optical path 
which turns out to be 1.56 m. 
3. We correct the geometrical PMT to source distance for the actual op- 
tical path. 

When doing this, the data scale as shown in Fig. [7} The fitting of the rescaled 
data set leads to an attenuation length of 216 ± 12 cm for our NdLS at [Nd] = 
6.5 g/1. 

A test with a longer cell and a full optical modeling of the light propagation 
is planned. 

To address the issue of NdLS stability, which is crucial among the others, 
we have kept the cell filled with NdLS along ~1 year and then performed a 
new set of measurements scanning as usual the cell along its main axis by the 
137 Cs source. In Fig. |8] we plot the quantity Ly® evaluated in 2009 spectra 
as a function of the source distance from each PMT; the only operation done 
on the NdLS prior to the measurements was nitrogen sparging to remove 
the oxygen that unavoidably entered in the cell, as the cell sealing is not 
guaranteed along such a long time. It is relevant to add that the laboratory 
where the cell was stored along the year was not air conditioned, therefore the 
temperature reached quite high values (~35°C along the summer months). 
The fitted uncorrected, for actual optical cell path and light propagation 
effects, A value is found to be 87±3 cm to be compared with 74±3 cm fitted 
from the May 2008 data set; therefore we can state that the scintillator is 
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Figure 6: Distribution of £f/f for mea- 
surement taken with the NdLS and 
normalizing data to corresponding PC 
measurements (see text for details). 



Data 2008 - lm cell 
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Figure 7: Distribution of L^E for measurements 
taken with the NdLS and using the actual cell 
optical path from timing data. 



stable along the observation period. 

To determine the light yield, we compute the quantity: 



PMTl ■ Qpmt2 oc LY ■ e A (2) 

where Qi is the charge of an event, LY is the light yield (p.e./MeV), L the 
length of the cell and A effective attenuation length for the 1-m cell. Then 
we plot it as a function of the source position taking as reference PMT2 as 
shown in Fig. [9] In particular, we are interested in the ratio: 

VQ^Si ■ 111 = LY NdLS ( ) 

V VPMTl VPA/T2 Lb 

The ratio in Eq ^ allows to determine the LY in the NdLS relative to that 
of the undoped LS, that is by how much the Nd reduces the scintillation 
light. This is shown in Fig. |9j where it can be seen that the NdLS reduces by 
~25% the light yield of the pure LS. This value is consistent with the value 
measured in the small vials at the time of production. 

4. Radiopurity and purification methods 

The intrinsic radiopurity of Nd has been measured with HP-Ge detec- 
tors at the Gran Sasso underground facility [9j. On a sample of 50.2 g of 
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Figure 8: Distribution of L^E for both PMTs against the relative distance to PMT2 (see 
text for detailes). 2009 data. 



commercial NCI2O3 interactions of photons produced by internal radioactive 
contaminants have been recorded for 37.3 days. In Table [T] we report the 
results of the measurement. 



228 Ra 


228 Th 


226 Ra 


234 Th 


234p a m 


235JJ 


40 K 


137 Cs 


ppb 


ppb 


ppb 


ppb 


ppb 


ppb 


mBq/k 


mBq/kg 


6.6±1.8 


8.5±1.1 


0.39±0.23 


<22 


<9.5 


<9.1 


<86 


9.7±2.7 



Table 1: Radiopurity measurement of commercial Nd 2 03 used for the NdLS production 
by means of HP-Ge. Upper limits are at 90% C.L. 



The most important source of background we are concerned about for 
applications of the NdLS is 232 Th. From data in Table [T] we conclude that 
in the commercial Nd 2 03 a contamination at the level of 8 ppb in 232 Th is 
present. This turns to about 30 mBq/kg. A second set of measurements 
were carried out by ICP-MS on different packages of the same production 
lot. Samples were diluted in a 10 ml HNO3 solution at 2.5%. In Table [2] we 
report results. 

Both methods applied to determine the intrinsic Th contamination in 
Nd2C>3 have shown similar results. In order to reduce the Th contamination 
we have developed and tested two purification methods. In a first case we 
have used a liquid-liquid purification. In particular, at room temperature we 
perform a liquid-liquid/ extraction using NdCl3 in Hydrochloric acid solution 
and an equal volume of TOPO diluted in toluene. Mixing lasted for about 
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Comparison PCLS vs NdLS 
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Figure 9: Distribution of L^E from the spectra built for ^Qpmi ■ Qpmti against the 
relative distance to PMT2 (see text for detailes). 



Nd sample 


104 mg 


237 mg 


0.067 mg 


232 Th [ppb] 


5±1 


7±2 


6±1 


238 U [ppb] 


< 0.2 


<0.3 


<0.3 



Table 2: Radiopurity measurement of commercial Nd2C>3 used for the NdLS production 
from ICP-MS. Upper limits are at 68% C.L. 



2 minutes. Using 0.2M TOPO in toluene and starting from a sample spiked 
with 200 ppb of Th we have reached a final contamination of 2 ppb, that is 
a purification factor of ~ 100. 

As an alterative to the liquid-liquid extraction method a selective ion 
exchange column procedure using UTEVA resin has been attempted. The 
active group in UTEVA resin is the dipentyl-pentyl -phosphonate that se- 
lectively extracts thorium and uranium from a solution containing an high 
(3M) concentration of nitric acid. A 15% Nd(NC>3)3 solution in 3M nitric 
acid artificially contaminated with Th is passed through a column filled with 
the UTEVA resin that retains Th and U. Immediately after, the column is 
washed with 3M nitric acid to recover most part of the Neodimium. The 
Neodimium loss is less than 1% and the purification factor, calculated mea- 
suring the Th concentration both in the un-purified and purified solution by 
mean ICP-MS, is about 200. The column can be used without significant 
loss of efficiency to purify about 15g of Nd for lg of resin. The column can 
be regenerated for some (less then 10) times using strong HC1 (6M) to elute 
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the thorium and the uranium from the column, then washing the column 
with high purity water and conditioning again the resin with 3M nitric acid. 
The advantages of this procedure are the higher purification factor that can 
be obtained and the fact that the method is well suited for a large-scale 
continuous operation. The disadvantages are the high cost of the resin and 
the fact that the Neodimium is recovered in a strongly acidic solution that is 
not usable for the Nd carboxylate synthesis. For this reason a further step is 
needed to neutralize the nitric acid, precipitate the Nd as hydroxide with the 
aid of urea, purify the precipitate from the nitric ion and finally to convert 
the hydroxide to chloride. This step complicate further the procedure so that 
only a marginal advantage remains from the significant (two-fold) increase 
of the purification factor. 

5. Neutrinoless double beta decay with a 12 tons Nd-loaded scin- 
tillator 

The neutrinoless double beta decay (Oz/2/3), (A, Z) — > (A, Z+2)+2e~ , will 
reveal the Majorana nature of neutrinos. For a given candidate nucleus (A,Z) 
the quantity to be measured is the decay rate: T^ 1 = G^M^rrii^ where Gi 
is a determinable phase-space factor, Mi the nuclear matrix element (NME) 
and m 2 pp the so-called effective Majorana neutrino mass. The quantity to 
be measured contains two unknowns. Therefore, at least two measurements 
are needed to determine tt&q. Extra measurements are needed to test the 
nuclear model. At present a claim of evidence comes from the search in 
76 Ge pU]. This claim turns into a 90% confidence interval for rrvig equal to 
[0.17,0.45] meV [UJ. Therefore, 100 meV can be taken as a reference value 
for coming experiments to probe the present finding. 

The idea to search for neutrinoless double beta decay with a Nd-loaded 
liquid scintillator is not new [12] . In this large mass of Nd is used in 

1 kton of liquid scintillator and this allows to search for such a rare decay 
process in spite of the poor energy resolution when comparing with other 
detectors [131 E3]- In the following we propose a similar measurement but 
using a much smaller detector. In particular, we aim to merge the Borex- 
ino technology [TJ)] , which has shown the feasibility to purify 1 kton of liquid 
scintillator at the level of 10~ 17 g/g in 238 U and 232 Th, with the results estab- 
lished in our measurements and considering the performances of new PMTs 
with high quantum efficiency (QE) and low radioactivity (< 100 mBq) [I6]. 
Taking into account the results obtained and reported in the previous Sec- 
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tions, we have studied the possibility to search for the Nd zero-is mode double 
beta decay with a 10-ton scale target mass of NdLS based on PC. In partic- 
ular, we can assume the same bulk radiopurity reached by Borexino because 
we are using the same solvent and fluor. Therefore, our basic assumptions 
are a PC-based organic liquid scintillator with 1.5 g/1 of PPO and 6.5 g/1 
of Nd. As described above this NdLS will have a light yield reduced by a 
25% with respect to pure PC-based scintillator and an attenuation length of 
~2 m. We make use of 12 tons of NdLS in a configuration similar to the 
one used in Borexino [15]. The NdLS is contained in a nylon vessel 1.5 m in 
radius inside a 4 m radius and 1 cm thick stainless steel sphere filled with 
distilled PC. As many as possible 3-inch low radioactivity and high quantum 
efficiency PMTs are installed on the inner surface of the sphere so that a 
80% of coverage is reached. Light concentrators can be used to reduce the 
total number of PMTs needed to get a large light collection [17\. A geomet- 
rical amplification factor on the order of 4 or larger can be assumed. Using 
4 as many as 2000 PMTs will be needed. The total activity at the inner 
nylon vessel due to PMTs is estimated to be ~ 2mBq, using 100 mBq/PMT 
and Apc = 27.2 cm for the absorption length of 2.6 MeV photons in PC. The 
stainless steel sphere will produce about 0.4 mBq using a steel contamina- 
tion of ~2 mBq/kg [18]. So, a maximum activity of about 2.4 mBq will be 
compared with the rate of the two-z/ decay mode. For such a detector we 
can determine the light yield, LY, as: LY = G\ x G c x QE x r\ x Q Nd x N y , 
where C7a=0.62 accounts for the reduction of light collection efficiency due 
to the finite attenuation lenght, C7 C =0.80 is the geometrical light collection, 
<5-E=0.4CQ r/=0.90 the PMT collection efficiency on the first dynode and 
iV 7 ~10 4 photons/MeV for a PC-based liquid scintillator. With the above 
assumptions we obtain at a first approximation about 1340 p.e./MeV and a 
resolution of 3.5% FWHM at (5/3,8=3.367 MeV. For Nd zero-v decay mode 
the half-life, T?Y 2 , in years can be written [19]: = ^4 — ml, where m v 

I 1/2 m e 

is in units of meV and the electron mass, m e , in units of eV. The quantity 
T?f 2 can be used to determine the expected rate for the zero-z/ decay mode. 



lr The quantum efficiency depends on the wavelength. Using data from |16j we have 
determined that the average value of QE is found at 395nm. However, the number of 
photoelectrons registered depend upon the QE, the emission spectrum and the trasmit- 
tance. As it is shown in Fig. [2] the NdLS has low attenuation between 370nm and 420nm, 
when this is taken into account the average sits at 391nm, where QE^0.40. Therefore, 
our assumption is well justified. 
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For the present study we do not consider uncertainties on predictions for the 
transition matrix related to the decay because this is beyond the purpose of 
this paper. 

In our target mass of 12 tons an intrinsic bulk contamination of Th at 
the level of lCT 5 ppt will produce about 6 cpy from 208 T1 (Qp=5 MeV). Only 
0.1 cpy will fall within la of Qpp. At a concentration of 6.5 g/1 about 92 kg 
of Nd are needed. For an effective Majorana neutrino mass of 300(100) meV 
we expect ~28(3) cpy due to the zero-z/ decay mode. The two-z/ decay mode 
gives about ~50 mBq using Ty 2 = 9.2 x 10 18 years. Therefore, this decay 
mode will give by far the largest counting rate below 3 MeV. Another im- 
portant source of background are 8 B solar neutrinos. These neutrinos will 
produce ~22 cpy and only about 0.3 are expected within la of Qpp. In 
Table [3] we summarize the expected rates for signal and backgrounds for 
771^=300(100) meV. In these two cases natural abundance (5.6%) is consid- 
ered. For the purpose of our study we have assumed a Th contamination of 
1(0.1) ppt in Nd and considered a 90% reduction due to 212 Bi- 212 Po tagging 
and statistical subtraction. In this scenario the Th contamination in the Nd 
will give about 10(1) cpd within la of Qpp. 





M Nd 




FWHM 


m v 


S 0v 


S 2 » 


8 B^ 


208 rjj 


s/Vb 


tons 


kg 


kg 


% 


meV 


cpy 


cpy 


cpy 


cpy 




12.3 


92 


5.1 


3.5 


300 


19.2 


0.1 


0.3 


10.6 


5.8 


12.3 


92 


5.1 


3.5 


100 


2.1 


0.1 


0.3 


1.1 


1.8 



Table 3: The case of a 12 tons detector loaded with 6.5 g/1 of Nd at natural abundance. 
All rates are given within la of Qpp- 

The case reported in Table [3] allows to measure the zero-v decay mode at 
the level of 17%(40%) in 3(5) years for m y =300(100) meV. 

Much better results can be obtained with Nd enriched in 150 Nd. An 
enrichment factor of 60% would give S' ^=23 cpy with 55 kg of 150 Nd within 
la and for m u =100 meV. Such a possibility could be within experimental and 
technological reach due to the limited 150 Nd mass needed [2D]. Moreover, for 
a small target mass the background due to solar neutrinos is negligible and 
not as important as it could be in the case of a larger detector such as SNO+. 
In conclusion with a target mass of 12 tons of LS, a 60% Nd enrichment and 
lppt of Th in Nd a measurement at the level of 15% (30%) seems feasible in 
3(5) years of data taking for m„=100(50) meV. In Tableland in Fig. 10 we 
summarize our results for the case of the enriched detector. 
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M L5 


M Nd 


Miso/vd 


FWHM 


m v 




S 2v 


8 B^ 


208rjJ 


s/Vb 


tons 


kg 


kg 


% 


meV 


cpy 


cpy 


cpy 


cpy 




12.3 


92 


55 


3.5 


100 


23 


1.0 


0.3 


10.6 


6.6 


12.3 


92 


55 


3.5 


50 


5.7 


1.0 


0.3 


10.6 


1.7 



Table 4: The case of a 12 tons detector loaded with 6.5 g/1 of Nd at 60% enrichment and 
lppt Th in Nd. All rates are given within la of Qpp- 
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Figure 10: Expected spectrum for the case of an NdLS enriched at 60% for 771^=100 meV. 



6. Conclusions 

Neutrinoless double beta decay is one of the most important research 
goals in Astroparticle Physics at present. In this paper we have shown how 
this search can be carried out by means of a highly doped liquid organic 
scintillator with a target mass of about 10 tons and using 150 Nd. We have 
reported results of measurements carried out to determine the optical proper- 
ties, time stability and photon yield for a 6.5 g/1 Nd doped liquid scintillator 
based on pseudocumene. In particular, with our experimental results we have 
shown that a target mass of about 12 tons with about 60 kg of enriched Nd 
could reach a sensitivity in the range of 50-100 meV in the effective neutrino 
mass. We have also reported about purification methods studied to reduce 
the intrinsic 232 Th contamination in commercial Nd. We believe our results 
could pave the way for further studies with the aim to better characterize 
the feasibility of the technique presented for neutrinoless double beta decay 
search. 
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